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SUMMARY 

Computer  codes  for  modeling  electric  fields  with¬ 
in  the  cavity  of  externally  sustained  electric  dis¬ 
charge  lasers  have  been  developed.  These  codes  use 
semidi rect/marchi ng  methods  in  body-fitted  coordinates 
and  are  fast  enough  to  make  interactive  design  practi¬ 
cal  . 

Introduction 

The  objective  of  the  computational  effort  des¬ 
cribed  herein  was  to  develop  computer  codes  for  rap¬ 
idly  and  accurately  modeling  the  electric  fields 
within  the  cavity  of  pulsed  electric  lasers  such  as 
the  VIPER  laser  (a  program  at  the  Air  Force  Weapons 
Laboratory).  These  codes  should  be  fast  enough  to 
make  the  interactive  design  process  practical.  The 
designer  should  be  able  to  perturb  the  laser  oper¬ 
ating  parameters  and/or  the  electrode  geometry,  and 
quickly  obtain  new  solutions,  graphically  presented. 

In  the  design  of  electron  beam  lasers,  it  is 
desirable  to  have  a  nearly  uniform  energy  deposition 
throughout  the  cavity.  This  energy  deposition  is 
governed  by  the  solution  of  the  non-linear  elliptic 
equation  for  electric  potential  <f>,  given  by 

V  •  aV4>  =  0  (1) 

where  the  conductivity  a  is  a  non-linear  function  of 
the  electric  field  E  =  .  The  solution  of  this 

equation  for  a  reasonable  grid  resolution  in  two 
dimensions  is  a  time-consuming  effort  using  conven¬ 
tional  methods.  Since  breakdown  depends  on  the  field 
near  the  electrode  surfaces,  accurate  representation 
of  the  surfaces  is  essential  for  realistic  modeling. 

Solution  of  the  Potential  Equation 
in  Cartesian  Coordinates  by  Semi  direct  Methods 

The  approach  used  in  the  current  work  is  to  solve 
this  equation  using  semidi  rect  non-linear  solution 
techniques  previously  developed  for  fluid  dynamics 
applications  by  the  present  author.  (See,  e.g..  Refs. 
1-3.)  The  equations  are  first  linearized  about  some 
initial  guess  for  the  solution.  The  linear  equations 
so  generated  are  generally  variable  coefficient,  non- 
separable,  linear  elliptic  equations.  These  are 
solved  using  marching  methods  for  linear  equations 
which  are  the  only  practical  way  (outside  of  brute 
force  direct  Gaussian  elimination)  for  solving  such 
non-separable  elliptic  equations  directly  rather  than 
iteratively.  A  sequence  of  linear  solutions  is  then 
used  in  a  quasi-Picard  iteration  to  solve  the  non¬ 
linearity.  Typically,  the  solution  of  the  non-linear 
equations  is  converged  wi thin  10-15  iterations. 


The  ionization  S  of  the  external  electron  beam 
gun  is  modeled  empirically,  following  Ref.  9,  by  the 
following  equation. 

-  atan  (2) 

where  y  =  E/2  *  V .  The  electron  beam  has  a  voltage  V 
and  a  width  2a  located  at  x  =  0  between  y  =  -a  and 
y  =  +a.  From  the  same  reference,  the  non-linear  con¬ 
ductivity  a  is  given  by 

a  =  C-E°-45S0,5,  C  =  0(1)  (3) 

Initially,  the  laser  cavity  was  modeled  in  a 
cartesian  coordinate  system  with  straight  electrodes. 

In  this  system,  we  can  obtain  completely  converged 
non-linear  solutions  in  approximately  5  seconds  on  a 
31x31  grid  using  a  time-shared  CDC  6600  computer. 

This  is  two  orders  of  magnitude  faster  than  the  com¬ 
puter  time  necessary  to  solve  the  problem  using  a 
triangular  finite  element  system. 

An  example  of  the  solution  obtained  with  this 
method  in  cartesian  coordinates  is  presented  in  Fig¬ 
ures  1  and  2.  Contour  plots  are  presented  for  the 
electric  field  E  and  the  energy  deposition  rate  =  oE2. 

Solution  of  a  Body-Fitted  Coordinate 
Transformation  by  Semi  direct  Methods 

The  next  step  in  this  research  was  to  obtain  solu¬ 
tions  of  the  non-linear  potential  equation  in  realistic 
cavity  geometries  with  non-rectangular  cathode  and 
anode  shapes.  It  is  sometimes  possible  to  obtain  such 
solutions  using  a  cartesian  coordinate  system  with  par¬ 
tial  cells  near  the  irregular  boundaries.  However,  it 
is  well-known  that  such  a  treatment  often  leads  to 
numerical  instabilities  and/or  numeri cal  inaccuracies. 
We  have  chosen  to  use  the  approach  of  using  a  numeri¬ 
cally  generated  body- fitted  coordinate  system,  follow¬ 
ing  Thompson,  et  al .  10  >u. 

In  the  body- fitted  coordinate  system,  the  physical 
domain  (x,y)  is  transformed  to  a  regular,  rectangular 
coordinate  system  in  coordinates  (e,n).  In  order  to 
insure  that  no  coordinate  lines  cross  in  the  trans¬ 
formed  (c,n)  plane,  the  coordinate  transformation  is 
generated  by  the  solution  of  a  coupled  pair  of  non¬ 
linear  elliptic  equations  for  x  and  y  given  by 
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S  =  exp(-yx)  atan 
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where 


a  5  x^2  +  y^2  (4c) 

B  =  x?xn  +  y5yn  (4d) 

Y  =  x?2  +  y^2  (4e) 

These  equations  are  solved  for  x(5,n)  and  y(£,n), 
with  the  boundary  values  on  x  and  y  being  specified 
by  the  boundary  shape  and  the  desired  mesh  density 
at  the  boundaries.  Second-order  accurate  finite 
difference  equations  are  used  throughout. 

Note  that  all  calculations,  including  the  calcu¬ 
lation  for  the  generation  of  the  coordinates,  are 
performed  in  the  (?,n)  plane  rather  than  the  physical 
(x,y)  plane.  Thompson,  et  al .  10  >n  use  a  point  SOR 
iteraton  scheme  to  solve  the  coordinate  system  trans¬ 
formation,  but  this  proves  to  be  inefficient  for  high 
mesh  resolution.  We  again  use  a  semi  direct  method 
to  solve  the  coordinate  system  transformation.  This 
system  is  well-suited  to  the  semi  direct  method  be¬ 
cause  the  two  non-linear  equations  for  x(6,n)  and 
y(S,n)  couple  only  at  interior  points.  (This  con¬ 
trasts  to  the  boundary  coupling  of  the  Navier-Stokes 
equations,  which  slows  convergence.  See  Refs.  1-3.) 

The  speed  of  convergence  of  these  non-linear 
equations  for  the  coordinate  systems  depends  on  the 
accuracy  of  the  initial  guess.  For  a  mild  problem 
in  which  the  electrode  shapes  are  a  slight  distortion 
from  rectangular,  an  adequate  initial  guess  is  ob¬ 
tained  by  simple  linear  interpolation  between  boundary 
values  of  x  and  y  in  the  (s,n)  plane.  In  these  cases, 
the  coordinate  system  transformation  is  solved  typi¬ 
cally  in  less  than  10  non-linear  quasi-Picard  itera¬ 
tions.  An  example  of  the  coordinates  so  generated  is 
shown  in  Figure  3.  For  more  severe  problems,  ade¬ 
quate  initial  guesses  are  difficult  to  achieve.  (For 
example,  even  for  a  simple  U-shaped  cavity,  linear 
interpolation  gives  a  folded  coordinate  system  with 
a  negative  Jacobian.)  For  these  cases,  it  will  be 
necessary  to  build  up  to  the  final  geometry  by  way 
of  a  sequence  of  intermediate  geometries,  using  non¬ 
linear  continuation  methods.  This  aspect  of  the  work 
is  still  under  development. 

Solution  of  the  Potential  Equation 
in  Body-Fitted  Coordinates 

The  transformed  non-linear  potential  equation  is 
now  to  be  solved  in  the  (?,n)  plane.  However,  the 
form  of  the  equation  changes  drastically,  because  the 
coordinate  transformation  used  is  non-conformal .  In 
particular,  mixed  or  cross  derivatives  of  the  form 
324>/358n  are  generated  where  none  existed  in  the 
physical  plane.  This  introduces  no  inaccuracy  in 
the  solution,  but  does  require  that  the  solution 
method  used  to  solve  the  equations  be  able  to  treat 
a  general  9-point  operator.  Again,  the  marching 
methods  for  elliptic  equations  are  capable  of  handling 
this  generality,  although  at  some  penalty  in  computer 
ti  me . 

For  simple  geometries  like  that  shown  in  Figure  3, 
in  either  planar  or  cylindrical  geometries,  the  non¬ 
linear  solution  on  a  25x25  grid  is  obtained  in  5  to 
15  seconds  on  a  CDC  6600. 

Interactive  Graphics 

An  important  goal  of  the  present  work  is  to  en¬ 
able  the  researcher  to  interactively  design  electrode 


shapes  and  laser  operating  parameters  so  as  to  achieve 
a  near  uniform  energy  deposition  the  laser  cavity. 

This  will  be  possible  because  of  the  rapid  solution 
time  of  the  semidirect  methods,  coupled  within  inter¬ 
active  software  currently  under  development. 

The  researcher  in  laser  physics  will  use  a  Tek¬ 
tronix  4662  plotter  to  digitize  the  electrode  and 
cavity  shapes  and  interactively  generate  body-fitted 
coordinate  systems,  then  interactively  generate  solu¬ 
tions,  perturb  the  parameters  of  the  problem,  and 
finally  obtain  acceptable  solutions  and  contour 
plots  of  the  solution  functionals. 

Future  Work 

In  the  immediate  future,  we  intend  to  develop  the 
three  areas  described  below. 

(1)  Coordinate  system  control  via  the  non- 
homogeneous  terms10’11  in  (4). 

(2)  Coordinate  solutions  for  severe  geometries 
via  non-linear  continuation. 

(3)  Interactive  design  capability  via  digitizing 
and  graphics  displays. 

In  the  longer  range,  we  intend  to  pursue  the 
following  areas. 

(4)  Improved  solution  for  non-linear  conduc¬ 
tivity,  possibly  using  a  Monte  Carlo 
solution. 

(5)  Improved  far- field  boundary  conditions, 
possibly  using  asymptotic  solutions. 

(6)  Three-dimensional  solutions. 

(7)  Fourth-order  accurate  solutions. 

(8)  One-dimensional  sheath  solution  of  the 
charge  separation  equations,  possibly  with 
Townsend  coefficients  and  diffusion,  cou¬ 
pled  to  the  multidimensional  electric  field 
codes . 

(9)  Automation  of  the  design  process  to  achieve 
a  near  uniform  energy  deposition. 

(10)  Time  dependent  solutions. 
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Figure  3.  Example  of  solution  of  body- fitted 
coordinates  for  a  mild  geometry. 
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Figure  1.  Example  of  electric  field  solution 
in  rectangular  coordinates. 


Figure  2.  Energy  deposition  =  aE2  for  the  electric 
field  solution  in  rectangular  coordinates. 
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